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Abstract. We show that the traditional criterion for a simplex to belong to the Delaunay 
triangulation of a point set is equivalent to a criterion which is a priori weaker. The 
argument is quite general; as well as the classical Euclidean case, it applies to hyperbolic 
and hemispherical geometries and to Edelsbrunner's weighted Delaunay triangulation. In 
^) \ spherical geometry, we establish a similar theorem under a genericity condition. The weak 

fNJ ' definition finds natural application in the problem of approximating a point-cloud data set 

with a simplical complex. 
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I. Strong and weak witnesses 

Let A C W 1 be a set of points, not necessarily finite. The Voronoi diagram of A, de- 
noted Vor(A), is the decomposition of M. n into Voronoi cells {V a : a G ^4} defined as follows: 

V a = {x G K n : \x - a\ < \x - b\ for all b G A} 



The dual of the Voronoi diagram is the Delaunay triangulation Del(A), an abstract simplicial 
complex which contains the p-simplex [ ] with vertices ao, Oi, . . . , a p G A if and only 

if V ao PI V ai fl • • • fl V ap ^ 0, for all p > 0. When A is a finite set of points in general position, 
there are no cells with p > n, and Del(A) is geometrically realised as a triangulation of the 
convex hull of A. In this paper our interest is in the abstract simplicial complex, which is 
defined under all circumstances. An example of the geometrically-realisable case is given in 
Figure [TJ the edges of the Delaunay triangulation of the set A = {a, b, c, d, e, /} are shown 
q . in black, and the Voronoi diagram is in grey. 

In order to state our theorem, we reformulate the definition of Del(A) in terms of witnesses. 

O ■ 

Definition. Let a = [ciq&i . . . a p ] be a p-simplex with vertices in A. We say that x G W 1 is 
a weak witness for a with respect to A if \x — a\ < \x — b\ whenever a G {ao, a>i, ■ ■ ■ , a p} and 
rS ! b G A \ {a ,ai, . . . ,a p }. 
c5 

Definition. Let a = [a Q ai . . . a p ] be a p-simplex with vertices in A. We say that x G M. n is 
a strong witness for a with respect to A if x is a weak witness for a with respect to A and 
additionally there is equality \x — a \ = \x — oi| = ••• = |x — a p |. 

Examples are given in Figure 121 In the left panel, a; is a strong witness for triangle [adf] . 
In the center panel, y is a weak witness for edge [ac] and indeed for triangle [abc\. On the 
right, z is a weak witness for [bed]. The circumcenter of [feed] lies close to a, so [bed] does not 
have a strong witness. 

When the context is unambiguous, we may occasionally drop the phrase "with respect 
to A" . It follows from the definition of Voronoi cells that a G Del(A) if and only if a has a 
strong witness. 
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Figure 1. The Voronoi diagram and Delaunay triangulation for a 6 point set. 






Figure 2. Strong and weak witnesses. 

Theorem 1. Let i C K" and let a = [a ai . . . a p ] be a simplex with vertices in A. Then a 
has a strong witness if and only if every subsimplex r < a has a weak witness. 

A consequence of this theorem is that the Delaunay triangulation can be defined using 
weak witnesses only. Note that one direction is trivial: if x is a strong witness for a then it 
is a strong (and hence weak) witness for every subsimplex r < a. 

Example. If abed is a convex quadrilateral in the plane, then each of the vertices [a], [b], [c], 
[d], each of the boundary edges [ab], [be], [cd], [da] and each of the triangles [abc], [bed], [cda], 
[dab] has a weak witness with respect to {a, b, c, d}. If abed is not cyclic, then exactly one of 
the diagonal edges [ac], [bd] has a weak witness. For instance, in Figure edge [ac] has a 
weak witness but edge [bd] does not. The theorem predicts that the two triangles containing 
that edge have strong witnesses; and indeed the Delaunay triangulation consists of those two 
triangles joined along a common edge. 



The special case of an edge a = [a ai] is well-known, and appeared in |MS94j to justify a 
construction of a topology-approximating graph. When their construction is generalised to 
topology-approximating simplicial complexes CdS03j, the corresponding justification is pro- 
vided by the full theorem as it appears here. The terminology of weak and strong witnesses 
was introduced in the latter paper. 



2. Proof of Theorem □ 

The proof of Theorem ^ is an inductive generalisation of the proof in [MS94 for the case 
of an edge. We strengthen the statement of the theorem, in order to facilitate the induction. 



Theorem 2. Let A cW 1 and let a = [a ai 



be a simplex with vertices in A. For every 



nonempty I C {ao, cl\, ■ ■ ■ , a p }, suppose there exists a weak witness xj for the subsimplex 
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Figure 3. Finding yy. the two possibilities in the case a\ £ J. 



aj < a with vertex set I. Then the convex hull of the points xj contains a strong witness 
for a. 

Proof. The proof proceeds by induction on p. The base case p = is trivial. Now suppose 
p > 1. Consider the subspace H = {x £ IR n : \x — ao| = \x — a\\} and the (p — l)-simplex 
t = \a\a2 ■ ■ ■ dp]. The following statements are easily seen to be equivalent: 

x £ H is a strong witness for a with respect to A 
x £ H is a strong witness for r with respect to A 
x £ H is a strong witness for r with respect to A' = A \ {ao} 

For every nonempty J C {ai, . . . ,a p }, we will find a weak witness yj E H with respect 
to A' for the subsimplex tj < r having vertex set J. By the inductive hypothesis, this 
implies the existence of a strong witness x for r with respect to A', in the convex hull of the 
points yi, and hence in H. The weak witnesses yj are constructed as convex combinations 
of the points xi, so x itself belongs to the convex hull of the points x\. This will complete 
the inductive step. It remains to locate the weak witnesses y,j. 

The first case is a\ £ J. Write J = {a\} U K where K C {02, . . . ,a p }, and consider 
weak witnesses x k, %\k and x ik for the simplices with vertex sets {a } U K, {a{\U K and 
{a ,ai} U K respectively. Suppose |x ia' — a i| — \ x oik ~ a o\- Since \xik — ai | < \xik — a>o\ 
it follows by continuity that there is some point y in the closed line segment £ = [x ik, Xik] 
at which \y — a\\ = \y — ao\. Moreover, for every a £ J and b £ A' \ J, the convex inequality 
\x — a\ < \x — b\ is valid at the endpoints of i and hence throughout I. It follows that yj = y 
is the required weak witness for tj. If instead \xoik — ai | < \xq\k — clq\, there is a similar 
argument with the roles of ao and a± interchanged, giving yj in [xoik,xok] H H. The two 
possibilities are illustrated in Figure El 

The second case is a% £" J, so 7^ J C {a 2 , . . . ,a p }. Consider weak witnesses xj, x j 
and xu for the simplices with vertex sets J, {a } U J and {ai} U J respectively. Suppose 
\%j — flol > — a i\- Since |x j — a \ < \x j — ai\, there is some point y in £ = [xj,x j] 
at which \y — ao| = \y — Moreover, for every a £ J and 6 £ A' \ J, the inequality 
|aj — a| < I a; — b\ is valid throughout I. Thus the required weak witness is yj = y. If instead 
\xj — ao| < |xj — ai|, we can interchange the roles of ao and a\ to find yj in [xj, ly] D H. 

This completes the proof of Theorem |21 and hence of Theorem ^ □ 

3. Variations 

The proof of Theorem |21 is quite general. Here is a set of sufficient requirements, in the 
form of axioms. 
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(1) Let R be a set with a notion of convexity; so every unordered pair x,y G R determines 
a subset [x, y] C R of intermediate elements. 

(2) Let A be a set. For every ordered pair (a, 6) in A, there is a convex subset R a b C -R, 
interpreted as "the set of points whose distance from a is no greater than their 
distance from b" . Convex here means closed under taking intermediate elements. 

(3) Define H ab = R ab fl R ba . If x G R ab and y G R ba then we require that [x, y] fl H ab ^ 0. 
This rule stands in for the intermediate value theorem. 

(4) For all a, b, c G A, we require that H ab fl R ac = H ab fl R bc and if a ft fl _R ca = H ab fl i? c ;,. 
In other words, a and f> can be treated alike on the subspace of points equidistant 
from them. 

One can define weak and strong witnesses in terms of the sets R ab , and formulate analogous 
versions of Theorems Q and 121 The proof goes through unchanged. Here are three important 
special cases. 

• Let R = TC n , hyperbolic n-space, and define R ab = {i£ TC n : d(x, a) < d(x, b)} using 
the hyperbolic metric. 

• Let R = S 1 ™, the open hemisphere in Euclidean n + 1-space, defined by 

SI = {x G W 1 : x\ + x\ + • • • + x 2 n = 1, x x > 0} 

and set R ab = {iG S™ : d(x, a) < d(x, b)} using the geodesic metric. 

• Let R = M. n and let a i— > w a be a real-valued function on A. Define R ab = {x G M" - : 
|x — a| 2 + w a < \x — b\ 2 + w b }. 

In each case, [x, y] is the set of points on the unique geodesic arc from x to y. The remaining 
axioms are easily verified; in particular the sets R ab are convex. It follows that weak witnesses 
can be used to define the Delaunay triangulation in hyperbolic space and the Euclidean 
hemisphere; as well as Edelsbrunner's weighted Delaunay triangulations |Ede95j in Euclidean 



space. In the latter case, it is not necessarily true that there is a weak witness for every 
0-dimensional simplex [a], so one must not fail to check this. In the other cases, a itself is 
always a witness for [a] . 

Remark. A stronger version of Axiom 4 is transitivity: R ab fl R bc C R ac . This holds in all 
our examples. 

4. The Euclidean sphere 

In this section and the next, we study the case of the Euclidean sphere S n , with weak 
and strong witnesses defined in terms of the intrinsic geodesic metric. The direct analogue 
of Theorem is false, so we must weaken the claim. We begin by reformulating Theorem ^ 
itself. 

Theorem 3. Let A C lR n and let a = [aoai • • • a p ] be a simplex with vertices in A. Suppose 
that for every subsimplex r < a there exists a closed Euclidean ball B T; such that B T contains 
the vertices of r and the interior of B T contains no other points of A. Then there exists a 
closed Euclidean ball B, such that the boundary of B contains the vertices of a and the 
interior of B is disjoint from A. 



Proof. The existence of B T is equivalent to the existence of a weak witness x T for r: take x T 
to be the center of B T , or conversely take B T to be the smallest ball centered on x T which 
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contains all the vertices of r. Similarly, the existence of B is equivalent to the existence of 
a strong witness for a. In this way Theorem El is equivalent to Theorem^ □ 

For the sphere S n we have the following result. 

Theorem 4. Let A C S n and let a = [aoai . . . a p ] be a simplex with vertices in A. Suppose 
that for every subsimplex r < a there exists a closed metric ball B r , such that B r contains 
the vertices of r and the interior of B T contains no other points of A. Suppose further that 
the union of the balls {B T : t < a} is not equal to the whole sphere S n . Then there exists a 
closed metric ball B, such that the boundary of B contains the vertices of a and the interior 
of B is disjoint from A. 

Proof. Let s G S n \ \J{B T : r < a}. Regard s as a vector in « n+1 and let V C R n+1 
be a hyperplane perpendicular to s. It is well known that the stereographic projection 
S n \ {s} — > V maps closed metric balls in S n which are disjoint from s onto closed Euclidean 
balls in V; and that this is a bijective correspondence. Under this correspondence, TheoremH] 
is equivalent to Theorem El □ 

Example. TheoremEJis false without the extra condition. Consider the poles n,s G S 2 . Let 
ai, a 2 , . . . , a p belong to a circle of latitude C, and let a be any point lying north of C. Every 
simplex r < a = [a ai . . . a p ] then has a weak witness: if Oq G t then n is a weak witness 
for r, and if Oq G" t then s is a weak (and indeed a strong) witness for r. The corresponding 
disks are the northern and southern polar caps bounded by C; these cover S 2 between them. 
If p > 3 then a has no strong witness. 

The proof of Theorem El breaks down because of the failure of convexity on S n . There 
is no unique shortest geodesic connecting a pair of antipodal points on S n . More seriously 
the O-sphere S° is disconnected, so there is not even one geodesic on which to find the 
intermediate point required by Axiom 3. The inductive argument works by reducing the 
problem to subspheres of successively lower dimension, so the S° case is unavoidable when 
p > n. On the other hand, for subspheres S with k > 1, the failure of convexity is non- 
generic, since it is only a problem when there is an antipodal pair of witnesses at some stage 
in the proof. We will develop this idea in Section 

Corollary 5. Suppose A = {a ,ai, . . . ,a n+ i} C S n is not contained in an (n — l)-sphere 
of any radius. For every J C A suppose there is a closed metric ball Bj containing J and 
whose interior is disjoint from A \ J. Then the balls Bj cover the sphere S n . □ 

Example. For < k < n, identify W n = R k © R n ~ k with the tangent space of S n C R n+1 
at the north pole. For 1 ^> e > e' > 0, let {oq, ai, . . . , a^} be the vertices of a regular 
/c-simplex on the sphere of radius e in M. k , and let {bo, b±, . . . , 6 n -fc} be the vertices of a 
regular (n — &)-simplex on the sphere of radius e' in IR n , projected onto a neighbourhood 
of the north pole in S n . If e' is sufficiently close to e, it can be checked that every simplex 
r < a = [a . . . Ofc&o • • • frn-fc] has a small witness disk B T contained in a neighbourhood of the 
north pole, except for r a = [a ai . . .a k }. Since a does not have a strong witness, it follows 
that every witness disk for the /c-simplex r a is necessarily large, and in particular contains 
the entire southern hemisphere. 

5. The Euclidean sphere, continued 



We will prove the following theorem. 
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Theorem 6. For a generic finite subset A C S n , and p < n, a simplex a = [aofli . . .a p ] 
with vertices in A has a strong witness if and only if every subsimplex r < a has a weak 
witness. Here "generic" means that, for all k < n, no set of k + 3 points in A is contained 
in a k-sphere (of any radius). 

Thus if A C S n is generic we can use weak witnesses to define its Delaunay triangulation 
in S n correctly as far as the n-skeleton. On the other hand, the true Delaunay triangulation 
is at most n-dimensional if A is generic: the existence of an (n + l)-cell would imply that 
its n + 2 vertices lay on a common (n — l)-sphere. The conclusion is that the Delaunay 
triangulation in S n of a generic set A is equal to the n-skeleton of the complex defined using 
weak witnesses. Since a randomly-chosen A C S n is generic with probability 1, this is a 
satisfactory result. 

Definition. Let A C S n (not necessarily finite), and let a ] be a p-simplex 

with vertices in A. We say that x G S n is a robust witness for a with respect to A if x has 
a neighbourhood in which every point is a weak witness for a with respect to A. 

The next proposition may clarify the situation. 

Proposition 7. For n > 1 , let A C S n be a finite set and let a = [a Oi . . . a p ] be a p-simplex 
with vertices in A. Then x is a robust witness for o if and only if there is strict inequality 
d(x, a) < d(x, b) whenever a G {a , a±, . . . , a p } and b G A \ {a , a±, . . . , a p }. 

Proof. Assume first that x is a robust witness, and o G {a , Oi, . . . , a p } and b G A \ 
{a , ai, . . . , a p }. Suppose d(x,a) = d(x,b) < tt. (The possibility d(x,a) = d(x,b) = n 
is ruled out since x has only one antipode.) Let y ^ x, b be any point on the minimal 
geodesic [x,b]. Then 

d(y, b) = d(x, b) — d(x, y) 

but 

d(y, a) < d(x, a) — d(x, y) 

since y is not on the minimal geodesic [x,a]. Thus d(y,b) < d(y,a) for points y which can 
be arbitrarily close to x. Since this would contradict the robustness of x it must be that 
d(x, a) < d(x, b) in all cases. 

In the converse direction, note that the inequalities d(x, a) < d(x, b) are open conditions 
and there are finitely many of them. If all of the appropriate inequalities are satisfied at x 
(making x a weak witness), then they are also satisfied on an open neighbourhood of x. Thus 
x is a robust witness. □ 

Theorem 8. Let A C S n (not necessarily finite), let p < n, and let a ] be a 

p-simplex with vertices in A. If every subsimplex t < a has a robust witness, then a has a 
strong witness. 

Proof. For k — 0, 1, . . . ,p, let o~k — [afca^+i . . . a p ], let A). — A \ {ai : i < k} and let be 
the great sphere defined by the equations: 

d(x, a ) = d(x, ai) = • • • = d(x, a^) 

Note that dim(H},) > n — k. We will show inductively that every subsimplex r < has 
a robust witness in H k with respect to A k . The case k = is our hypothesis, and the case 
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Figure 4. [left] The n-plane V determined by the sphere C; [right] the sep- 
arating hyperplane W determines a suitable direction for perturbation. 

k = p implies the conclusion, since a witness x G H p (robust or otherwise) for [a p ] with 
respect to A v is automatically a strong witness for a. 

Suppose the assertion is proved for some k < p. For each r < cr^ we have a robust 
witness x T G H]~ with respect to A^. We may assume that no two of the x T are antipodal; since 
dim(iTfc) > 1 we can perturb the points x T to make this true. We now apply the argument in 
the proof of Theorem El to find weak witnesses y T G H k+1 for r < cr fc+1 with respect to A k+ i. 
If an inequality d(x, a) < d(x, b) is satisfied on neighbourhoods of points xo, x\ which are not 
antipodal, then it is satisfied on a neighbourhood of the minimal geodesic interval [xo,Xx]. 
Thus the new witnesses y T are themselves robust. This completes the inductive step. □ 

If A is generic we can manufacture robust witnesses from weak witnesses. 

Lemma 9. Suppose A C S n is a finite set satisfying the condition that, for all k < n, no 
set of k + 3 points in A is contained in a k-sphere. Let a be a simplex with vertices in A. If 
o has a weak witness, then a has a robust witness. 

Proof. Let x be a weak witness for a. If x is not robust, then there are points ai, ci2, . . . , a& 
which are vertices of a, and b\, 62, • • • , G A which are not, such that 

d(x, ax) — ■ ■ ■ — d(x, au) = d(x, b\) = ■ ■ ■ = d(x, h), 

so all these points lie on a (n — l)-sphere C centered on x. The other vertices of a may 
be assumed to lie inside C, and the remaining points of A outside C. By the condition, 
k + l < n + 2. 

Regarding S n as a subspace of IR n+1 , consider the n-plane V defined by V H S n = C. 
This can naturally be identified with the tangent space of S n at x. We will construct a 
vector w G V, and show that, for a small perturbation x' = x + ew + 0(e 2 ), there are strict 
inequalities d(x',ai) < d(x',bj) for all i,j. When e is small enough, the relevant inequalities 
for the remaining points of A remain unchanged, so x' will be the required robust witness. 

Inside V, we claim that the convex hull of {ai, . . . , a^, 61, . . . , b{\ is an embedded (k+l — 1)- 
simplex. If this were false, then the k + I points would lie on a (k + / — 2)-plane in M. n and 
hence on a (k + I — 3)-sphere in S n , contradicting the hypothesis. It follows that the convex 
hulls of {ai, a,2, • • ■ , au] and {61, 62, • • • , h} are disjoint, so they can be strictly separated by 
an (n — l)-plane W . Let w be the unit normal vector of W < V, pointing towards the 
half-space containing the points aj. See Figure EJ 

We now determine the effect of the perturbation x' = x + ew + 0(e 2 ) along S n , on the 
distances. It is simplest to compute the derivative of the squared Euclidean distance function 
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/ a :6H \x + ew — a| 2 , namely /„(0) = —2(w,a) (since (w,x) = 0). Since w was chosen as the 
normal to a separating hyperplane, we deduce the inequality (0) — f' b . (0) = 2(w, bj—cii) < 
for all i, j. Thus for small e we have |x' — | 2 < \x' — bj\ 2 and hence c?(x', eij) < 6 3 -) since 
there is a monotonic relationship between the two distance functions. This completes the 
proof. □ 

Proof of Theorem® We have the following sequence of implications: 

a has a strong witness =>- every r < a has a weak witness 
[by Lemma 01 generically] =^ every r < a has a robust witness 
[by Theorem |Hl since p < n] =>• a has a strong witness 
This proves the theorem. □ 

6. Application 

It is a standard problem in computational geometry to recover the topology of a manifold 
M C R n from a discrete sample of points A C M. Recovery is often taken to mean 
the construction of a simplical complex with vertices in A, which is homeomorphic to the 
unknown space M. This problem has a rich literature. 

When M is regarded as known, it is natural to consider is the restricted Delaunay trian- 
gulation Del(A, M), which is dual to the partitioning of M by Voronoi cells. Unlike the full 
Delaunay complex Del(A), which is always contractible, Del(yl, M) carries information about 
the topology of M and may even be homeomorphic to it, if A is sampled sufficiently finely. 
According to the traditional definition, a simplex a with vertices in A belongs to Del(v4, M) 
if there is a strong witness for a on M. 

When M is not known, one may attempt to approximate Del(A,M). This approach, 
introduced in jMS94j for approximation by graphs, is extended to higher dimensional com- 
plexes in |CdS03j . The idea is to select a subsample L C A to serve as the vertex set, and to 
let the remaining points stand in for the unknown manifold M. Unfortunately, it is useless 
to define Del(L, A) by mimicking the traditional definition of Del(L,M), since strong wit- 
nesses in A exist with probablility 0. The solution is to replace definitions involving strong 
witnesses with definitions involving weak witnesses, which do exist with nonzero probability. 

The goal of this paper has been to establish the credentials of weak witness definitions, 
proving equivalence in the archetypal cases of Delaunay triangulation in Euclidean, hyper- 
bolic and spherical geometries, and paying attention to some of the pitfalls. All this may 
be regarded as groundwork for the following question: under what conditions is the weak- 
witness complex Del w (L, A) homeomorphic (or homotopy equivalent) to M? 
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